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ABSTRACT 

As shown by Parker and Raval, quantum field theory in curved spacetime 
gives a possible mechanism for explaining the observed recent acceleration of 
the universe. This mechanism, which differs in its dynamics from quintessence 
models, causes the universe to make a transition to an accelerating expansion in 
which the scalar curvature, R, of spacetime remains constant. This transition 
occurs despite the fact that we set the renormalized cosmological constant to 
zero. We show that this model agrees very well with the current observed type- 
la supernova (SNe-Ia) data. There are no free parameters in this fit, as the 
relevant observables are determined independently by means of the current cosmic 
microwave background radiation (CMBR) data. We also give the predicted curves 
for number count tests and for the ratio, w(z), of the dark energy pressure to 
its density, as well as for dw(z)/dz versus w(z). These curves differ significantly 
from those obtained from a cosmological constant, and will be tested by planned 
future observations. 

Subject headings: cosmic microwave background — cosmological parameters - 
cosmology: observations - - cosmology: theory - - gravitation — supernovae: 
general 

1. Introduction 

Observational evidence appears increasingly strong that the expansion of the universe 
is undergoing acceleration that started at a redshift z of order 1 (Riess et al. 1998, 2001; 
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Perlmutter et al. 1998, 1999). Observations of scores of type-la supernovae (SNe-Ia) out to 
z of about 1.7 support this view (Riess et al. 2001), and even glimpse the earlier decelerating 
stage of the expansion. It is fair to say that one of the most important questions in physics 
is: what causes this acceleration? 

One of the more obvious possible answers is that we are observing the effects of a small 
positive cosmological constant A (Krauss & Turner 1995; Ostriker & Steinhardt 1995; Do- 
delson, Gates, & Turner 1996; Colberg et al. 2000). Another, less obvious possibility is that 
there is a quintessence field responsible for the acceleration of the universe (Caldwell, Dave, 
& Steinhardt 1998; Zlatev, Wang, & Steinhardt 1999; Dodelson, Kaplinghat, & Stewart 2000; 
Armendariz-Picon, Mukhanov, & Steinhardt 2001). Quintessence fields are scalar fields with 
potential energy functions that produce an acceleration of the universe when the gravita- 
tional and classical scalar (quintessence) field equations are solved. More recently, Parker & 
Raval (1999a,b,c, 2000, 2001) showed that a quantized free scalar field of very small mass 
in its vacuum state may accelerate the universe. Their model differs from any quintessence 
model in that the scalar field is free, thus interacting only with the gravitational field. The 
nontrivial dynamics of this model arises from well-defined finite quantum corrections to the 
action that appear only in curved spacetime. This was the first model to present a realiza- 
tion of dark energy with ratio of pressure to energy density taking values more negative than 
— 1. Other models having this property have subsequently been proposed (Caldwell 2002; 
Melchiorri et al. 2002). 

The physics of the Parker-Raval model is based on quantum field theory in curved 
spacetime. The renormalized (i.e., observed) cosmological constant A is set to zero. Several 
mechanisms have been proposed that tend to drive the value of A to zero (Dolgov 1983; Ford 
1987, 2002; Tsamis & Woodard 1998a,b; Abramo, Tsamis, & Woodard 1999), but these 
mechanisms play no role in this model. The energy-momentum tensor of the quantized field 
in its vacuum state is determined by calculating an effective action (Schwinger 1951; DeWitt 
1965; Jackiw 1974) in a general curved spacetime. The spacetime is unquantized, and is 
itself determined self-consistently from the Einstein gravitational field equations involving 
the vacuum expectation value of the energy-momentum tensor, as well as the classical energy 
momentum tensor of matter (including cold dark matter) and radiation. In solving the 
Einstein equations, the symmetries of the FRW spacetime are imposed, but the spacetime is 
not otherwise taken as fixed, and the initial value constraints of general relativity are satisfied. 
The acceleration is the result of including in the effective action a non-perturbative term 
involving the scalar curvature of the spacetime (Parker & Toms 1985a,b,c; Jack & Parker 
1985). The minimal effective action that includes this non-perturbative effect and gives the 
correct trace anomaly of the energy-momentum tensor was used by Parker and Raval. In 
applying this effective action to the recent expansion of the universe, terms involving more 
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than two derivatives of the metric were neglected. In this approximation, a solution was 
proposed in which the universe undergoes a rapid transition from a standard FRW universe 
dominated by cold dark matter to one containing significant contributions of vacuum energy 
and pressure. The proposal is that this negative vacuum pressure is responsible for the 
observed acceleration of the universe. The reaction back of this negative vacuum pressure 
on the expansion of the universe is such as to cause a rapid transition to an expansion of the 
universe in which the scalar curvature remains constant. The vacuum pressure and energy 
density are determined by a single parameter related to the mass of the scalar field. The 
transition to constant scalar curvature is the result of a rapid growth in the magnitudes of 
the vacuum pressure and energy density that occurs in this theory when the scalar curvature 
approaches a particular value, of the order of the square of the mass of the particle associated 
with the scalar field. The Einstein equations cause a reaction back on the metric such as to 
prevent further increase in the magnitudes of the vacuum pressure and energy density. (This 
effect is analogous to Lenz's law in electro magnetism.) 

The essential cosmological features of this model may be described quite simply. For 
times earlier than a time tj (corresponding to z ~ 1), the universe undergoes the stages of the 
standard model, including early inflation, and radiation domination followed by domination 
by cold dark matter. During the latter stage, at time tj the vacuum energy and (negative) 
pressure of the free scalar quantized field increase rapidly in magnitude (from a cosmological 
point of view). The effect of this vacuum energy and pressure is to cause the scalar curvature 
R of the spacetime to become constant at a value Rj. The spacetime line element is that of 
an FRW universe: 

ds 2 = -dt 2 + a(t) 2 [(l - kr 2 Y l dr 2 + r 2 d6 2 + r 2 sin 2 6d(j) 2 ] , (1) 

where k = ±1 or indicates the spatial curvature. By joining at tj the matter dominated 
scale parameter a(t) and its first and second derivatives to the solution for a(t) in a constant 
R universe, one uniquely determines the scale parameter a(t) for times after tj. This model 
is known as the vacuum cold dark matter (VCDM) model of Parker and Raval, or as the 
vacuum metamorphosis model (Parker & Raval 1999c) to emphasize the existence of a rapid 
transition in the vacuum energy density and pressure. 

The constant value, Rj, of the scalar curvature is a function of a single new parameter, 
m, related to the mass of the free scalar field. Therefore, the function a(t) for t > tj is fully 
determined by fh. The values of m and tj can be expressed in terms of observables, namely, 
the present Hubble constant H , the densities Q m0 and Q r0 of the matter and radiation, 
respectively, relative to the closure density, and the curvature parameter Qko = — ^/(-^o°o)- 
(Here a = a(t ) is the present value of the cosmological scale parameter.) These observables 
have been determined with reasonably good precision by various measurements that are 
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independent of the SNe-Ia (Krauss 2000; Freedman et al. 2001; Hu et al. 2001; Huterer & 
Turner 2001; Turner 2001; Wang, Tegmark, & Zaldarriaga 2002). Therefore, the value of m 
is known to within narrow bounds, independently of the SNe-Ia observations. 

The power spectrum of the CMBR depends largely on physical processes occurring long 
before tj. The behavior of a(t) in the VCDM model does not significantly differ from that 
of the standard model until after tj. Therefore, the predicted power spectrum in the VCDM 
model differs only slightly from that of the standard model. We calculate the predicted 
power spectrum of the CMBR in the VCDM model (as described below), and find the range 
of values of the above observables that give a good fit to the CMBR observations. 

From this range of observables, the corresponding range of the parameter m follows. 
Therefore, the prediction of the VCDM model for the magnitude versus redshift curve of 
the SNe-Ia is completely determined, with no adjustable parameters. We plot the predicted 
curves (obtained from this range of fh) for the distance modulus A(m — M) of the SNe-Ia 
as a function of z. Comparison with the observed data points, as summarized by Riess et 
al. (2001), shows that a significant subset of predicted curves fit the SNe-Ia data very well, 
passing within the narrow error bars of each of the binned data points, as well as of the 
single data point at z ~ 1.7. 

We also give the curves predicted by the VCDM model for number counts of cosmological 
objects as a function of z and for the ratio, w(z), of the vacuum pressure to vacuum energy 
density, as well as for dw/dz versus w (parametrized by z). The predictions of the VCDM 
model differ significantly from those of the ACDM model. Accurate measurements of these 
quantities out to z of about 2 would be very telling. 

The model we consider here is the simplest of a class of models in which a transition 
occurs around a finite value of z within the range of possible observation. This general class 
has been studied from a phenomenological point of view by Bassett et al. (2002a,b). They 
find that the CMBR, large scale structure, and supernova data tend to favor a late-time 
transition over the standard ACDM model. Although they considered only dark-energy 
equations of state, w, with w > —1, their phenomenological analysis can be generalized 
to include the present VCDM model. In addition, the VCDM model is readily generalized 
to include a nonzero vacuum expectation value of the low mass scalar field, which could 
bring w into the range greater than —1. In the present paper, we are taking the simplest 
of the possible VCDM (or vacuum metamorphosis) models, so as to introduce no arbitrary 
parameters into the fit to the supernova data. 
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2. 



How Observables Determine m 



In this section, we explain how the value of m is obtained from Hq, Qko, ^mOj and 
Q r o in the VCDM model. (Here, Q m o = fi c dmo + ^&o, where fi c dmo and Qbo are the present 
densities of cold dark matter and baryons, respectively, relative to the closure density.) The 
relation between rh and these present observables follows from the Einstein equations, the 
previously described constancy of the scalar curvature R = Rj for t > tj, and continuity of 
a{t) and its first and second derivatives at time tj. For our present purposes, we may define 
the parameter rh in terms of Rj, namely, by the relation rh 2 = Rj. (At a microscopic level, 
rh is proportional to the mass of the free quantized scalar field.) 

The trace of the Einstein equations at time tj is 



where p m j is the energy density of the non-relativistic matter present at time tj. The density 
and pressure of the dark energy, p v and p v , respectively, will be taken to be zero at tj. 
We make this assumption in order to avoid introducing a second parameter in addition to 
rh. It should be noted for future reference that dropping this assumption will introduce 
another parameter that would affect mainly the behavior of the predicted SNe-Ia curve near 
the transition time tj from matter-dominated to constant-scalar-curvature universe. In the 
present paper, we do not relax this assumption of zero p V j and p V j because we find good 
agreement of the one parameter VCDM model with the current observational data. 

For all t > tj, the scalar curvature is taken to remain constant at the value Rj = rh 2 . 



m 2 = Rj = 8nGp. 



(2) 



Thus, 




(3) 



where dots represent time derivatives. Defining the variable x = a 



2 , this becomes 



-x + k — -m 2 x . 



(4) 



The first integral of this equation is 



-x 2 — —m 2 x 2 + kx = E , 



(5) 



where E is a constant. 



One of the Einstein equations at time tj is 




(6) 
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where Hj = a(tj)/a(tj), a,j = a(tj), and pj = p(i,) is the total energy density at time t = tj. 
The last equation can be rewritten as 



-x 2 + kxj = (8nG/3)pjX 2 j , 



(7) 



where subscript j refers to quantities at time tj. Comparing this with equation (5), one finds 
that 

lnG/3)p j - \-m 2 



E 



12 „. (8) 

Using pj = p m j + p r j, where p r j is the radiation energy density at time tj, and equation (2) 
to eliminate pj and p m j from the last expression for E, we find that 



E = (87rG/3)p rj x 2 + -m 2 x 2 . 



Thus, equation (5) gives the following conserved quantity: 



11 1 

-x 2 — —m 2 x 2 + kx = (8nG/3)p r jX 2 + -m 2 x 2 . 



12 



(9) 



(10) 



This is readily written in terms of a(t) and its derivatives. We have x/x — 2a/ a = 2H(t), 
and x 2 /x(t) 2 = [aj/a(t)] 4 . Hence, the radiation energy density satisfies 



p rj x 2 /x(t) 2 = p r (t) , 



and equation (10) is 

H(tf + k/a(t) 2 = 
Solving for [a J /a(t)] 4 , we obtain (for t > tj), 



irG/3)p r (t) + ^-m 2 [a 3 /a(t)f + ^-m 2 



K'/ a (^)] 4 = ^2 [ H ^f + k/^? ~ (8vrG/3)p r (t)] - \ 



(11) 



(12) 



(13) 



Returning to equation (2), we can now express m in terms of the present values of p m , 
p r , H, and k/a. We have 

rh 2 = Rj = 87iGp m0 (a /aj) 3 . (14) 



With equation (13), it follows that 

m 2 = 8nG Pm0 [Hi + k/a 2 - (87rG/3)p r0 ] - ^| 



-3/4 



(15) 
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Using the expression for the present critical density, p c0 = 3Hq/(8itG), equation (15) takes 
the dimensionless form, 



where mno = m/H . This is readily solved numerically for m 2 m . Alternatively, it can be put 
into the form of a fourth-order equation for m 2 m and solved analytically. Using the values 
O m0 = 0.34^^4 and H = 65+}q km s _1 Mpc -1 obtained in section 5 from the CMBR power 
spectrum (see fig. 2), as well as f2 r0 = 8.33 x 1CT 5 and il k0 = 0, we find m H o — ^-^Q+o'll 
and m = 4. x 10~ 33 eV. (The uncertainties refer to the 95% confidence level.) Note 
that the fitting of the CMBR power spectrum alone gives us a wide range for Q m o and 
H . In order to further restrict our results (thus being able to make stronger predictions), 
we will adopt the HST-Key-Project result H = 72 ± 8 km s' 1 Mpc -1 (Freedman et al. 
2001) as a constraint. This narrows the range of our cosmological parameters down to 
Q m0 = 0.34+°;^ and H = 65~^ km s" 1 Mpc" 1 , as can be seen from figure 2. Aiming at 
combining these two methods of determining the uncertainties of our results, hereafter we 
adopt the notation exemplified by fi m0 = 0.34^[ ) 1 4 6 ' ' 08 ' ) and H = 65^ 6 '^ km s _1 Mpc -1 , 
where the uncertainties appearing in parenthesis refer to the 95% confidence level, without 
and with the HST constraint, respectively. (The sign appearing in front of the parenthesis 
is common to both uncertainties.) Thus, returning to the parameters m m and m, we have: 




3/4 



(16) 



m H0 = 3.26 



-(0.62;0.08) 
'+0.10 



(17) 



and 



m = 4.52 



,-(1.76;0.18) 
+0.84 



x 10 



33 



eV . 



(18) 



Finally, the solution to equation (4) for x(t) = a(t) 2 is 




(19) 



It then follows from (l/2)x(t)/x(t) = H(t), that 




(20) 
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From equations (12), (14), and (16) one obtains the following expression in terms of the 
redshift, z = a /a — 1: 



H(z) 2 /H 2 



kO 



m 2 m /12)(l + zf + n k0 (l + zf + m 2 m /l2 . 



From equations (14) and (17) one finds the redshift Zj at time tj, 

Ko/(3^o)] 1/3 -i = i.i 9 ;(™- 19) . 



Z 3 



(21) 



(22) 



Moreover, from equation (21) we can obtain the redshift, z a , at which the expansion of the 
Universe starts to accelerate. In fact, by looking at the deceleration parameter, 



aa 



- = (l + z) [\n(H/H )]'-l, 



(23) 



where the prime sign stands for derivative with respect to z, we have that z a satisfies 

(l + z a )H'(z a ) = H(z a ) . (24) 

Thus, from equation (21) and the cosmological parameters mentioned above, we obtain for 
the spatially flat VCDM model 



Z a 



m 2 m /12 



1 - Q k0 -m 2 m /12 



1/4 



- 1 = 0.67 



-(0.59;0.15) 
+0.35 



(25) 



This value is similar to the one obtained using the spatially flat ACDM model with 
0.67: 



ZaACDM 



2Q 



AO 



1 — ^A0 



1/3 



0.60 . 



(26) 



3. The Dark Energy 

In the VCDM model, the dark energy is the energy of the vacuum, denoted by p v . This 
vacuum energy is not in the form of real particles, but may be thought of as energy associated 
with fluctuations (or virtual particles) of the quantized scalar field. Vacuum energy, p v , and 
pressure, p v , must be included as a source of gravitation in the Einstein equations. Thus, 
for t > tj, one has 

H{tf + k/a{t) 2 = (8ttG/3) \p r {t) + p m (t) + p v (t)} . (27) 



The vacuum energy and pressure remain essentially zero until the time tj when the value 
of the scalar curvature R has fallen to a value slightly greater than rh 2 . Then in a short 
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time (on a cosmological scale), the vacuum energy and pressure grow, and through their 
reaction back cause the scalar curvature to remain essentially constant at a value just above 
Rj = m 2 (Parker & Raval 1999a,b). Intuitively, this reaction back may be thought of as 
similar to what happens in electromagnetism when a bar magnet is pushed into a coil of wire. 
The current induced in the coil produces a magnetic field that opposes the motion of the 
bar magnet into the coil (Lenz's Law). Similarly, in the present case, the matter dominated 
expansion of the universe causes the scalar curvature to decrease. But as it approaches the 
critical value, m 2 , the quantum contributions to the energy-momentum tensor of the scalar 
field grow large in such a way as to oppose the decrease in R that is responsible for the 
growth in quantum contributions. The universe continues to expand, but in such a way as 
to keep R from decreasing further. 

Defining tj as the time at which p v and p v begin to grow significantly, we have to good 
approximation, equation (2). Evolving p mj forward in time, then gives 

8nGp rn (t) = m 2 K-/a(t)] 3 . (28) 

One then finds from equation (27) and equation (12) that the vacuum energy density evolves 
for t > tj as 

Pv{t) = 32^G t 1 " 4 [a ^ /a(t)]3 + 3 [°X*)1 4 } • ( 29 ) 
The conservation laws for the total energy density and pressure and for the energy densities 
and pressures of the radiation alone and of the matter alone, then imply that p v and the 
vacuum pressure, p v , also satisfy the conservation law. It follows that 

m 2 r , r , /.s-,4 



32nG 



{-l + [aX0] 4 } • (30) 



At late times, as a,j/a(t) approaches zero, one sees that the vacuum energy density and 
pressure approach those of a cosmological constant, A = m 2 /4. But at finite times, their 
time evolution differs from that of a cosmological constant. 

One immediately sees from equations (28) and (29) that for t > tj, 

— 2 

Pv{t)+ Pm {t) = ^{l + 3[a, J /a{t)f} . (31) 
Using equation (30), it now follows that 

p v (t) = (1/3) [p v (t) + p m (t)\ - m 2 /(24nG) . (32) 
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Since p m = 0, and p r = (l/3)p r , the total pressure, p, and energy density, p, satisfy the 
equation of state (Parker & Raval 2000) 

p(t) = (l/3)p(t) - m 2 /{2A7iG) . (33) 

From this equation of state and the conservation law, which can be written in the form, 
d(pa 3 ) I da + pd(a 3 ) / da = 0, we find that for t > tj, 

"' a »-(«-3^)<°>' 4+ 3^' < 34 > 

As the vacuum energy density is taken as zero at t — tj, we have pj = p m j + Prj = Pmoo^/a 3 + 
PrQa^/a^, and then, from equation (14), 

Pj = rh 2 /(8nG) + p r0 [m 2 /(87rGp m0 )] 4/3 . (35) 

Then, with equation (28) and (for t > tj) (a^/a) 4 = (ao/a) 4 (87rG ? p m o/m 2 ) 4 ^ 3 , equation (34) 
finally becomes 

p(a) = [pmo (3 W4)) V3 + Pro] (a /a) 4 + m 2 /(327tG) . (36) 

This expression will be used in the next section to calculate, among other things, the age of 
the Universe as predicted by the VCDM model. 

4. Age of the universe 

The values of tj and to are found by integration: 

ft fa(t) 

t= dt= daa- 1 H(a)- 1 . (37) 

Jo Jo 

From equation (27), 

pa(t) 

t= da a' 1 {-k/a 2 + (8irG/3)[p r (a)+ p m (a)+ p v (a)]}~ 1/2 , (38) 

Jo 

This integral is conveniently split in two at time tj, and expressed in terms of variable of 
integration y = a/ao'. 

fo-j/ao 

tj = (Ho)- 1 / d y (n k0 + Qmoy' 1 + Vroy- 2 ) ' , (39) 

Jo 
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and from equation (36), 

to - tj = {Ho)- 1 j dy {Q k0 + [Q m0 (3fi m0 /^ ) 1/3 + fi ro] V~ 2 + {m 2 m /12) y 2 }' 1 ' 2 . 

(40) 

Another interesting parameter to obtain is t a , the time when the expansion of the Universe 
starts to accelerate: 



to - t a = (#o) _1 / dy {Q k0 + [Q m0 (3ft m o/™D 1/3 + fi ro] V' 2 + {m 2 m /12) y 2 } 1/2 , 

J a a /a 

(41) 



where a a = a(t a ). 



Using the cosmological parameters obtained in section 5 by fitting the VCDM model 
to the CMBR power spectrum data, we have (see eqs. [22] and [25]) aj/a = (1 + Zj)- 1 = 
0.46^Q° g 4 '°' 04 ' ) and a a /a = (1 + £ a ) _1 = 0.60^q° 1 q 2 '°' 06 ' ) , and consequently, H tj = 
0.354^f - 009) , H t a = O.SSS^? - 013 ), and H t = O^; 2 ^. Thus, using the corre- 
spondent values of H , we finally obtain 

t = 5.33tS 8 T 22) Gyr , (42) 

t a = S-O-.-f"-' Gyr , (43) 

and 

to = 14.9 =F 0.8 Gyr . (44) 

Note that the uncertainty in t is independent of whether or not we adopt the HST-Key- 
Project result as a constraint. This is because the dependence of to on H and Q m o is 
such that the 95% confidence region shown in figure 2 happens to be stretched along lines 
of constant values of t . The same is true for the ACDM model, as discussed by Knox, 
Christensen, & Skordis (2001). 

We can compare the values presented in equations (43) and (44) with the respective 
ones given by the spatially flat ACDM model. Using il A0 = 0.67, the ACDM model gives 
H t a 0.536 and H t ~ 0.938. With the best-fit value of the Hubble constant obtained for 
our model, H = 65 km s _1 Mpc -1 (see section 5), this gives t a w 8.1 Gyr and t ~ 14.1 Gyr. 
Therefore, we see that the age attributed to the Universe by the VCDM model is larger than 
the age predicted by the ACDM model, for essentially the same values of f2 m0 and H . 



5. Fit to the CMBR power spectrum 



In view of recent CMBR observations (Pryke et al. 2002; Masi et al. 2002; Netterfield et 
al. 2002; Abroe et al. 2001), there is a need to reexamine the results obtained by Parker & 
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Raval (2001). In this section, we will obtain the cosmological parameters u c dmo = ^cdmo^ 2 , 
uj b0 = Q,b h 2 , and if = 100/i km s -1 Mpc -1 which give the best fit of the spatially flat VCDM 
model to the recent measurements of the CMBR power spectrum by Boomerang, maxima, 
and dasi (h is a dimensionless quantity defined by the latter expression). As we have seen, 
these parameters are the essential ingredients necessary to fix the parameter m. [The other 
necessary parameter, Q r o = 8.33 x 10~ 5 , is independently obtained from the CMBR mean 
temperature and the number of relic neutrino species (Peebles 1993).] 

In order to obtain the CMBR power spectrum fluctuations predicted by the VCDM 
model, with given values of u; c d m o, ^bo, and Hq, we use a slightly modified version of the 
cmbfast 4 computer code (Seljak & Zaldarriaga 1996; Zaldarriaga & Seljak 2000). The mod- 
ifications made in the code, described by Parker & Raval (2001), consist of adding the 
vacuum contributions, p v and p v , to the total energy density and pressure, respectively, for 
time t > tj, i.e., after the transition to dark energy dominated epoch. 

We set up the cmbfast code to generate a numerical grid in the 3 dimensional cosmo- 
logical parameter space (u> c dm0> ^bo, H ). We introduce prior information on the value of 
the present day Hubble constant, H , to be in the range from 45 to 80 with units of km 
Mpc -1 s _1 . Also, we set the value of the cosmological constant, A, to be zero. To per- 
form our numerical analysis consistent with these priors, we generated a class of VCDM 
models with the following cosmological parameters and resolutions (in the form of "initial 
value" : "final value":"step size"): u m = (0.005 : 0.030 : 0.001), u cdm0 = (0.05 : 0.31 : 0.01), 
and H = (45.0 : 80.0 : 1.0) km s" 1 Mpc -1 . We chose to vary these three parameters 
based on the fact that they determine directly m, which is the one free parameter of the 
VCDM model. All models generated use the Radical Compression Data Analysis Package 5 
(RadPack) (Bond, Jaffe, & Knox 2000) to compute a x 2 test statistic that compares the 
predicted CMBR spectrum to the experimental measurements of dasi (Pryke et al. 2002), 
Boomerang (Masi et al. 2002; Netterfield et al. 2002), and maxima (Abroe et al. 2001) at 
particular multipoles I. We look for minima of x 2 m the class of cosmologies specified above. 
The particular VCDM model described by the parameters which give the minimum of x 2 m 
our parameter space is called the best fit model. 

For the best fit VCDM model, we found oj m = 0.022, u cdm0 = 0.12, and H = 65.0 
km s _1 Mpc -1 . This best fit has x 2 = Xmin = 48.84, corresponding to a significance level 
aixmin) = g f(x 2 ,n) d( X 2 ) = 0.187, where f( X 2 ,n) = [2r(n/2)]- 1 (x 2 ^/ 2 ^-^ is 
the x 2 probability density function (pdf) for n degrees of freedom (the recent CMBR data 



4 http://physics.nyu.edu/matiasz/CMBFAST/cmbfast.html 
5 http: / /bubba.ucdavis.cdu/^knox/radpack.html 
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Fig. 1. — Plot of the best fit CMBR power spectrum for the spatially flat VCDM model (solid 
curve) and of the CMBR power spectrum for the spatially flat ACDM model with f2 A o = 0.67 
(dashed curve). The diamond points and error bars correspond to dasi, Boomerang, and 
maxima experimental data. 
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we consider are constituted of 41 experimental points). The three best fit cosmological 
parameters of the VCDM model are similar to the values obtained from the ACDM model 
fit to the CMBR power spectrum (Pryke et al. 2002; Abroe et al. 2001; Masi et al. 2002; 
Netterfield et al. 2002). 

To find the 95% confidence region of our parameter space, we compute the quantity 
Xg 5% where a(x95%) — 0.05, and then look for the subset of VCDM models with x 2 < Xg^% 
which lie in the parameter space given above. The results for the VCDM model are H = 
65^ Jo km s" 1 Mpc -1 , co c d m o = 0-12ioo3> an d ww = 0-022^qoo4- The uncertainties given above 
are obtained from the extrema of the 95% confidence region. The range of is consistent 
with Big Bang Nucleosynthesis (Bean, Hansen, & Melchiorri 2002), and the best fit value 
for H lies within the range given by the HST Key Project (Freedman et al. 2001). 

From the 95% confidence region in the parameter space (c^ c dmO) H ) we can obtain 
the 95% confidence region in the space (Q m o, H ), recalling that Q m0 = (u; c d m o + ^bo)/h 2 
and h = H / (100 km s" 1 Mpc" 1 ). This confidence region is presented in figure 2 (dark-gray 
region). The reason to concentrate in the parameter space (fl m o, H ) is that all quantities we 
are interested in (e.g., the age of the Universe, the luminosity distances of SNe-Ia, number 
count tests, the dark energy equation of state) depend directly on Cl m0 , not on cj cdm0 and u b0 
separately (except, obviously, the quantities directly related to the CMBR power spectrum). 
From figure 2 we see that Q m Q = 0.341°'^. This very wide range is just a consequence 
of the fact that the Hubble constant H is not tightly constrained by the CMBR power 
spectrum. However, adopting the HST-Key-Project result Hq = 72 ± 8 km s" 1 Mpc" 1 as a 
constraint (light-gray region in fig. 2), the acceptable range of Q m o can be narrowed down to 
O m0 = 0.34^^4, which then allows us to make stronger predictions about the cosmological 
quantities mentioned above. Using the even stronger prior H = 65 km s" 1 Mpc" 1 (the 
best-fit value) we obtain the "best-fit range" f2 m0 = 0.34 ± 0.06, which will be used in the 
next section to fit (with no free parameters) the luminosity distances of SNe-Ia. 

For the CMBR power spectrum predicted by the best fit flat VCDM model (see fig. 1), 
the multipole numbers / and power intensities // = 1(1 + \)Ci/{2n) of the first three peaks 
and two troughs are given in table 1. The uncertainties Al and A// correspond to the 95% 
confidence region of u cdmQ , u b0 , and H . Comparing the values given in table 1 with the 
correspondent results presented by Durrer, Novosyadlyj, & Apunevych (2001), we see that 
our values for the first three peaks and two troughs are well within the (la) ranges for these 
quantities that can be obtained, in a model-independent way, from the combined Boomerang, 
maxima, and dasi data. For instance, the last two columns of table 2 of Durrer, Novosyadlyj, 
& Apunevych (2001) give the multipole number and the power intensity of the first peak to 
be l Pl = 213+59 and L pi = 5041 ±ii96 /"K 2 ; respectively. 
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Fig. 2. — Plot of the 95% confidence level (dark-gray region) in the parameters Q m o and H 
obtained by fitting the CMBR power spectrum using the VCDM model. The extrema of 
this region, together with the best-fit values, define the ranges Q m0 = O.SA^fl and H = 
65;}q km s" 1 Mpc -1 . Adopting the HST-Key- Project result (H = 72 ± 8 km s" 1 Mpc" 1 , 
represented by the light-gray region), the range in Q m0 is narrowed down to Q m o = 0.34 + 
(dotted lines). Fixing H Q = 65 km s _1 Mpc -1 , we have f2 m o = 0.34 ± 0.06 (dashed lines) 
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The results of this section show that the VCDM model gives a reasonable fit to the 
CMBR power spectrum with values of a; c dm(b ^bo, and H that are consistent with current 
observations. The future of CMBR observations looks very promising with a mixture of 
ground based interferometers (dasi), airborne interferometers (maxima and Boomerang), and 
satellite experiments (Microwave Anisotropy Probe and the Planck satellite) that will further 
probe the CMBR anisotropies at higher and lower multipoles. 



6. No-parameter fit to the SNe-Ia data 

In the present section we compare the luminosity distance as a function of redshift 
predicted by the VCDM model to the measured values of the luminosity distances of SNe-Ia 
as summarized by Riess et al. (2001). Because all the relevant parameters of the model 
are determined by fitting the CMBR power spectrum (see sec. 5), this comparison is a 
no-parameter fit of the VCDM model to the SNe-Ia data. 

We start by computing the luminosity distance as a function of redshift, di(z). As a 
consequence of equation (1), the comoving coordinate distance r(z) of objects observed with 
redshift z satisfies 



dr dt dz 



VI - kr 2 a(t) a H(z) ' 

which leads to 



(45) 



J^aoHizT'dz' ,k = 

" 1 5 ^ < [a H n%) 1 sinh [H Qli 2 £ H^dz) , k + ' (46) 



where H(z) for the VCDM model is given by (see eq. [21]) 

H(z) = | [(1 - f4o - m 2 H0 /12) (1 + zf + n k0 (1 + zf + m 2 m /\2] 1/2 ,z< z 3 

H ° 1 [n r o (i + zf + n m0 (i + zf + n k0 (i + zf] 1/2 ,Z> Z 



(47) 



J 



The luminosity distance and the distance modulus are defined respectively as 

d L (z) =a (l + z)r(z) (48) 

and 

A ( m -M)(z) = 51og(|^, (49) 

where d L1 is the luminosity distance in the spatially flat VCDM model and d L2 is the lu- 
minosity distance in an arbitrary fiducial model used as normalization. We will set d^2 as 
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Table I. Multipole numbers I and associated power intensities I\ of the first few peaks 
and troughs of the CMBR spectrum predicted by the spatially flat VCDM model (see 
fig. 1). The uncertainties A/ and A J; are calculated from the x 2 test statistic. 
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Fig. 3. — Plot of the type-la supernovae distance modulus (normalized to a spatially open and 
empty cosmos) as a function of redshift z for the spatially flat VCDM (with fl m0 = 0.34±0.06) 
and ACDM (with Q A0 = 0.67) models. 
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the luminosity distance in an open and empty Universe [a(t) = t and k = — 1], which is the 
convention used by Riess et al. (2001). 

It is important to note that the expression for A (m — M) does not depend explicitly 
on the present value of the Hubble constant H . However, if we adopt the results of the 
previous section summarized in figure 2, then fixing different values of H in that figure 
lead to different 95%-confidence-level ranges for Q m o, which in turn give rise to different 
predictions for the SNe-Ia luminosity distances. In particular, using for H the best-fit value 
Hq = 65 km s _1 Mpc^ 1 gives us f2 m o = 0.34 ± 0.06 (see fig. 2). In figure 3 we plot the 
distance modulus as a function of redshift predicted by the VCDM model using this "best- 
fit range" for Q m o, as well as the observed distance moduli of SNe-Ia. It can be seen from 
figure 3 that fixing the values of f2 m0 and H that best fit the CMBR data also gives a 
very good no-parameter fit to the SNe-Ia data. Moreover, any value of f2 m0 in the "best-fit 
range" f2 m0 = 0.34 ± 0.06 gives a reasonably good fit to the SNe-Ia data, as shown by the 
dashed curves in figure 3. We also show in figure 3 the distance modulus predicted by the 
ACDM model with Q^o = 0.67. We see that even though the predictions of both models 
differ significantly in the range 0.5 < z < 1.5, current data are still not able to make a clear 
distinction between them. 

More numerous and accurate data on SNe-Ia luminosity distance are expected for the 
near future. The planned Supernova Acceleration Probe (SNAP) 6 , for instance, aims at 
cataloging up to 2,000 SNe-Ia per year in the redshift range 0.1 < z < 1.7. This improvement 
in our knowledge of the luminosity distances of SNe-Ia will provide a much stronger test of 
the VCDM model. 



7. Number Counts 

Counting galaxies or clusters of galaxies as a function of their redshift seems to be a 
very promising way to test different cosmological models (Huterer & Turner 2001; Podariu 
& Ratra 2001). The idea behind this procedure is that once we know, either by analytic 
calculations or by numerical simulations, the evolution of the comoving (i.e., coordinate) 
density of a given class of objects (e.g., galaxies or clusters of galaxies), counting the observed 
number of such objects, per unit solid angle as a function of their redshift, is equivalent to 
tracing back the area of the Universe at different stages that we can observe today. In other 
words, it is equivalent to determining our past light cone by constructing it from the area of 
these observed spherical sections of the Universe, parametrized by their redshift. Since this 



6 http://snap. lbl.gov/ 
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light cone is very sensitive to the underlying cosmological model, number counts provide a 
valuable tool for testing the mechanism which accounts for the accelerated expansion of the 
Universe. 

This kind of test was first performed using galaxies brighter than certain (apparent) 
magnitudes by Loh & Spillar (1986), with the simplified assumptions that the comoving 
density of galaxies is constant and that their luminosity function retains similar shape over 
the redshift range 0.15 < z < 0.85. Using the photometric redshift of 406 galaxies in that 
range, they were able to measure the ratio of the total energy density in the Universe to the 
critical density, obtaining f2o — 0.9^q 5 . However, the validity of Loh & Spillar's assumptions 
is still not clear due to the lack, to the present, of a complete theory of galaxy formation and 
evolution. In order to circumvent this problem, Newman & Davis (2000) then suggested that 
galaxies having the same circular velocity may be regarded as good candidates for number 
count tests, since the evolution of the comoving number density of dark halos having a given 
circular velocity can be calculated by a semi-analytic approach. Moreover, they claim, the 
comoving abundance of such objects at redshift z ~ 1 (relative to their present abundance) 
is very insensitive to the underlying cosmological model (under reasonable matter power 
spectrum assumptions). Other objects that one can count are clusters of galaxies (Bahcall 
& Fan 1998; Blanchard & Bartlett 1998; Viana & Liddle 1999; Haiman, Mohr, & Holder 
2001; Newman et al. 2002), which are simpler objects than galaxies, in the sense that their 
formation and evolution, and therefore their density, depend mostly on well-understood 
gravitational physics. 

Whatever class of objects one uses to perform the number count test, a key ingredient 
one needs to provide as an input, as stressed above, is the evolution of their comoving density, 

y/T=kFdN{z) 
HC{Z) = r 2 drdn ' (50) 
where dVl = sin OdOdcj) is the solid-angle element and dN(z) is the number of such objects, 
at the spatial section at redshift z, contained in the coordinate volume r 2 drdVt/ y/l — kr 2 . 
In order to find out the number of objects, per unit solid angle, with redshift between z 
and z + dz, we have to use the fact that the objects we are observing today with redshift 
z possess coordinate r which satisfies the past light-cone equations (45) and (46). Thus, by 
making use of these equations to eliminate the explicit radial dependence in equation (50), 
we get the number of observed objects with redshift between z and z + dz, per unit solid 
angle: 



dN j n c (z) (a H y 3 E(z)- 1 (f* E(z')- l dzf dz , k = 

dldn^ dz = j n c {z) (a H y 3 E(z)- 1 [fi~ 1/2 sinh (7^ 2 f?^')" 1 ^)] 2 dz ,k^0 



(51) 
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with E(z) = H(z)/H . 

In order to illustrate number counts predicted by the VCDM cosmological model, in 
figure 4 we plot dN (z) / dzdVt given by equation (51), with H(z) given by equation (47). Then, 
we apply the resulting formula to the spatially flat case and use the parameters obtained 
in section 5 adopting the HST-Key-Project constraint (see fig. 2), namely f2 m0 = 0.34^^4. 
Also, we follow Podariu & Ratra (2001) and Loh & Spillar (1986) in assuming, for simplicity, 
the constancy of the comoving density, n c (z) = noflg ( n o is the proper density at the present 
epoch). For sake of comparison, we also plot the spatially flat ACDM prediction, with 
^ao — 0.67 and the same assumption of constant density. Obviously, the predictions could 
be improved by dropping this latter assumption and taking into account the density evolution 
of the observed objects, as mentioned earlier. However, calculating such evolution is beyond 
the scope of the present paper, not to mention the fact that it is still not completely clear 
which class of objects we should choose. Moreover, once a more precise n c (z) is known, it 
is a simple task to take it into account since dN(z)/dzdfl is simply proportional to n c (z). 
From figure 4 we see that the VCDM model predicts more objects to be observed over 
redshifts z < 2 than the ACDM model. In fact, in a small redshift interval Az <C 1 around 
z ~ 1 the VCDM model predicts approximately 30% more objects than the ACDM model 
for approximately the same value of Q m0 (and the same value of n / Hq) . Note that this last 
conclusion should also hold for the counts of galaxies at fixed circular velocities suggested 
by Newman & Davis (2000), since their comoving density, even though not constant, is very 
insensitive to the underlying cosmological model at z ~ 1. We did not mention here the 
presence of selection effects, since they highly depend on the measurement procedure itself. 
Notwithstanding, these effects may also be included in the computation via an "effective" 
n c (z), which then should be viewed as the number of objects at the spatial section with 
redshift z, per comoving volume, satisfying the detectability conditions. 

Measurements of dN/ dzdQ will provide a valuable way to test the VCDM cosmological 
model and distinguish it from the ACDM model, when combined with CMBR anisotropy 
results. Such measurements will soon become available, as the DEEP (Deep Extragalactic 
Evolutionary Probe) Redshift Survey 7 expects to complete its measurements of the spectra 
of approximately 65, 000 galaxies in the redshift range 0.7 < z < 1.5 by the year 2004. 



7 http://deep. ucolick.org/ 
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Fig. 4. — Plot of the predicted number counts of objects per redshift interval and per solid 
angle, dN(z)/dzdQ (normalized by the present number of such objects in the volume Hq 3 ) as 
functions of their redshift z, for the spatially flat VCDM (with tt m0 = 0.34±g;2|) and ACDM 
(with Qaq = 0.67) models. In this plot, the comoving density is assumed to be constant. 
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8. Vacuum Equation of State 

The dark-energy equation of state p v = p v {Pv) in the VCDM cosmological model with 
k — ±1 or can be easely obtained, for t > tj, from equations (29) and (30): 



Pv = 3p v + 



— 2 

m 

8^G 



32nG ^ 3/4 
1 - ( 1 + — — Pv 



(52) 



Moreover, from the same pair of equations we also obtain the ratio w = p v /p v as a function 
of redshift: 



w(z) 



c 4 -i_ 

+- 1 

< C < 1 , (53) 



3(4 - 4C 3 + 1 

C 3 + C 2 + C + i 



3C 3 - c 2 - c - i ' 



where ( = aj/a = (1 + z)/(l + Zj). Note that equations (29), (30), (52), and (53) are the 
same as the respective ones presented by Parker & Raval (2001) in dealing with the spatially 
flat VCDM model. (Note, however, that the spatial curvature changes the value of Zj\ see 
eqs. [22] and [16].) 

In figure 5 we plot the redshift dependence of the ratio w, given by equation (53), for 
the spatially flat VCDM model using the cosmological parameters obtained in section 5 
adopting the HST-Key-Project constraint. The present value of this ratio, wo = w(z = 0), 
using Q m0 = 0.34^q° 1 4 6 '°' 08 ' ) , is w = — 1.28+o° 1 2 1 ' ' 08 ' ) . Note, from the expression for w(z), 
that w — > — oo as z — > Zj-, which is simply a consequence of the previously mentioned fact 
that the vacuum energy approaches a negligibly small value (more rapidly than the vacuum 
pressure) That no drastic consequence follows from the divergence of w is evident 

from figure 6, where we plot, using equations (33) and (36), the ratio w to t = p/p between 
the total pressure p and the total energy density p present in the Universe, as a function of 
redshift. For the sake of comparison, we also plot in the same figure the correspondent ratio 
Wfot given by the spatially flat ACDM model, 

(A) Pr + PA 



Pr + Pm + PA 

(lW3)(l + ^) 4 -^ A0 

O. (i + zy + (i - n r0 - n A0 )(i + zf + Ovo ' 



(54) 



with f^Ao = 0.67 and f2 r0 = 8.33 x 10~ 5 . In the VCDM model note that for times earlier than 
tj (i.e., z > Zj), during the matter dominated era, the total energy density, p ~ p m , and the 
total pressure, p = p r /3, lead to a negligible value of the ratio w to t- After tj (z < Zj), the 
negative pressure of the vacuum grows very rapidly in magnitude, becoming dominant and 
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Fig. 5. — Plot of the predicted ratio w = p v / p v as a function of redshift z, for the spatially 
flat VCDM model with Q m0 = 0.34l ^4 (the HST-Key-Project constraint was adopted). The 
present value of such ratio is wq = — 1.28^o!i2- Moreover, w — > — oo as z — > Zj- = (1. 19^47)- 
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determining a very sharp transition to the dark-energy dominated era. This is an important 
distinction between the VCDM (i.e., vacuum metamorphosis) model and the ACDM model, 
which presents a rather gradual transition (see fig. 6). 

In order to analyze not only the value of w but also its rate of change in redshift, we plot 
in figure 7, using f2 m0 = 0.34, the curve (w(z),w'(z)), where w'(z) = dw(z)/dz. The redshift 
z is used as the parameter of the curve. The present value w' = w'(z = 0) predicted by the 
spatially flat VCDM model is w' = -0.8 + ( q A ' '. Note that w -> -1 and w' -> as z — > — 1, 
which means that in the asymptotic future (assuming that nothing new will prevent the 
unbounded expansion of the Universe) the dark energy of the VCDM model behaves like an 
effective cosmological constant, with value given by A cff = m 2 /4 = 5.i;£ 2;0 ' 4) x 1(T 66 eV 2 . 
For comparison, in a ACDM model with H 65 km s" 1 Mpc -1 and f2 A o 0.67, one would 
have A = 3fi A0 # 2 « 3.8 x 10~ 66 eV 2 . 

The experimental determination of w(z), avoiding model-dependent assumptions, relies 
basically on measurements that, at least in principle, will determine H(z) with sufficient pre- 
cision to provide also a reliable determination of H'(z) = dH(z)/dz (Huterer & Turner 2001). 
To see this, let us consider the conservation equation satisfied by the total energy density p 
and the total pressure p, namely d(pa 3 ) + p d (a 3 ) = 0. Thus, with the (only) assumption 
that matter and radiation are separately conserved, we have that the energy density, px, 
and pressure, px, of dark energy (whatever it is) also satisfy the same conservation equation, 
which implies 

1 dp x _ (l + w x )d(a 3 ) _ 3 (1 + wx) ^ 



px dz a 3 dz (1 + z) ' 

where wx = Px/ Px and we have used 1 + z = a /a. Considering the general expression for 
the Hubble parameter as a function of redshift (which is obtained from Einstein's equation 
together with the assumption of separate conservation of matter and radiation), 

= n r0 (i + z f + n mQ (i + z f + n k0 (i + zf + n xo ^- (56) 

[with pxo being the present value of the dark-energy density and Qxo = 87rGpxo/(3i^o)L 
and using it and its redshift derivative to evaluate the left-hand-side of equation (55), we 
finally obtain the desired expression for wx(z): 

on (A (1 + ^) [dp x {z)/dz] 

wx{z = — ^ 1 

3 px{z) 



1 [2(1 + z)E'(z) - 3E(z)} E(z) - (1 + zf [ft r0 (l + zf - Q 



feO 



3 E{zf - n r0 (i + zf - n m0 (i + zf - n k0 (i + z 



(57) 
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Fig. 6. — Plot of the ratio w tot = p/p as a function of redshift for the spatially flat VCDM 
model with fi m0 = 0.34+^ and for the spatially flat ACDM model with fi A0 = 0.67. With 
these parameters, the present values of w to t and w[^t are w to t{z = 0) = —0.85^q^j and 
w[^(z = 0) = —0.67. Note that for the VCDM model the ratio w tot is negligible for times 
earlier than tj (i.e., z > Zj), during the matter dominated era. 
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Fig. 7. — Plot of the curve (w(z),w'(z)), parametrized by the redshift z, for the spatially 
flat VCDM model with Q m0 = 0.34. The solid-line curve corresponds to redshift z > 
while the dashed-line one corresponds to — 1 < z < 0. At the present epoch we have 
(wq,w' q ) = (—1.28,-0.8). Note that (w(z),w'(z)) — > (—1,0) as z — >• —1, which means that 
the dark energy of the VCDM model behaves like a cosmological constant in the asymptotic 
future. 
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where, again, E(z) = H(z)/H and E'(z) = dE(z)/dz. Thus, as stated above, wx(z) can 
be found from the determination of #(2) and H'(z). The quantity i^(^) can be determined 
from the direct observables d^z) and dN(z)/dzdfl, and the quantity n c (z) (Huterer & Turner 
2001). This is done by using equation (45) to express the derivative with respect to r in 
equation (50) in terms of a derivative with respect to z, and then using equation (48) to 
express r(z) in terms of di{z). This leads to the following expression for H(z): 



Thus, by considering measurements of luminosity distances and number counts, H(z) can 
be regarded as a directly observable quantity, which gives wx(z) through equation (57). In 
this sense, future data provided by the proposed satellite SNAP on supernovae luminosity 
distances (see sec. 6) and by the DEEP redshift survey on number counts (see sec. 7) may 
greatly improve our knowledge of the dark-energy equation of state. 



We have shown that the current observational data indicating that the expansion of the 
Universe is undergoing acceleration are quite consistent with the hypothesis that a transition 
to a constant-scalar-curvature stage of the expansion occurred at a redshift z ~ 1 in the 
spatially flat FRW universe having zero cosmo logical constant. This is the scenario proposed 
in the VCDM (or vacuum metamorphosis) model introduced by Parker and Raval. The late 
constancy of the scalar curvature at a value Rj = fh 2 is induced by quantum effects of a free 
scalar field of low mass in the curved cosmological background. The parameter fh, related 
to the mass of the field, is the only new relevant parameter introduced in this model, and 
can be expressed in terms of the present cosmological parameters Hq, f2 m0 , f2 r0 , and Qko (see 



Comparison of the CMBR-power-spectrum data with the flat-VCDM-model prediction, 
without or with the HST-Key-Project result as a constraint (see figs. 1 and 2, and table 1), 
gives the values of the cosmological parameters to be Hq = 65 + [q 6 ' 1 ' ) km s" 1 Mpc -1 and Q m o = 
0.34] 

(0^46,0.08)^ (R eca j] the definition of our notation in sec. 2: the uncertainties appearing 
in parenthesis refer to the 95% confidence level, without and with the HST constraint, 
respectively.) Such values lead to fh — 4.52+q 1 8 4 6,0 ' 18 - ) x 10~ 33 eV, and the best-fit values 
from the CMBR data give rise to a very good no-parameter fit to the SNe-Ia observational 
data. However, the SNe-Ia data are not accurate enough to draw a clear distinction between 
the VCDM and ACDM models. Other quantities of interest predicted by the VCDM model 
with the cosmological parameters mentioned above are the time and redshift at the transition 




(58) 



9. Conclusion 



eq. [16]). 
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between the matter-dominated and constant-scalar-curvature stages, (tj = 5.33^0.84 ' Gyr 
and Zj = 1.19^o°47 6 ' ' 19 ' ) ), the time and redshift when the accelerated expansion started (t a = 

g_ +(5.2;0.4) Qyr ^ = 67 -(O59;0.15)^ ^ ^ ^ q{ ^ Universe; ^ = 14 g Q.8 Gyr. 

Regarding future tests of the VCDM model, we have presented the prediction of number 
counts as a function of redshift, and compared it with the analogous ACDM prediction (see 
fig. 4). For approximately the same cosmological parameters, the VCDM model predicts 
nearly 30% more objects to be observed in a small redshift interval around z « 1 than 
the ACDM model. Data provided by the DEEP Redshift Survey in the near future will 
likely be able to distinguish these two models. Also, DEEP data combined with future 
measurements of SNe-Ia luminosity distances provided by the proposed SNAP satellite should 
greatly improve our knowledge of the dark energy equation of state, which bears the most 
distinct feature of the VCDM model: w < — 1 and w' < (see figs. 5 and 7). 

It should be noted that we have here considered the simplest form of the VCDM model, 
in which the transition to constant scalar curvature is continuous and effectively instanta- 
neous (see fig. 6). This form of the model makes definite predictions regarding the distance 
moduli of SNe-Ia and number counts. Thus, it is encouraging that it remains a viable model 
when confronted with the current observational data. Other natural parameters that may 
come into the VCDM model are the time interval over which the transition occurs, and the 
vacuum expectation value of the scalar field. A nonzero value of the transition time interval 
would mainly affect the predictions around z ~ 1, and a nonzero value of the vacuum ex- 
pectation value is likely to increase the ratio of pressure to density, w. Future observational 
data will determine if it is necessary to consider nonzero values for these parameters. 

This work was supported by NSF grant PHY-0071044 and Wisconsin Space Grant 
Consortium. The authors thank Koji Uryu for helpful comments and suggestions on the 
figures, and Alpan Raval for helpful discussions. 
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